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Abstract 

The hyperspherical harmonic basis is used to describe bound states in an j4-body system. The 
approach presented here is based on the representation of the potential energy in terms of hyper- 
spherical harmonic functions. Using this representation, the matrix elements between the basis 
elements are simple, and the potential energy is presented in a compact form, well suited for nu- 
merical implementation. The basis is neither symmetrized nor antisymmetrized, as required in 
the case of identical particles; however, after the diagonalization of the Hamiltonian matrix, the 
eigenvectors reflect the symmetries present in it, and the identification of the physical states is 
possible, as it will be shown in specific cases. We have in mind applications to atomic, molecular, 
and nuclear few-body systems in which symmetry breaking terms are present in the Hamiltonian; 
their inclusion is straightforward in the present method. As an example we solve the case of three 
and four particles interacting through a short-range central interaction and Coulomb potential. 

PACS numbers: 31.15.xj, 03.65. Ge, 36.40.-c, 21.45.-v 
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I. INTRODUCTION 



The Harmonic Hyperspherical (HH) method provides a systematic way to construct an 
expansion basis for a system of A particles. The N Jacobi vectors (N = A — 1) are trans- 
formed to the hyperradius p plus 3N — 1 hyperangular coordinates which are used to define 



' the hyperangular part of the 
1] and references therein). 



the HH functions. These functions are the eigenfunctions o 
Laplacian operator for a given number of particles (see Ref. 

Applications of the HH method to describe bound states of A = 3, 4 nuclei are well 
documented in the literature (for a recent review see Ref. |2j). In these applications the HH 
basis elements, extended to spin and isospin degrees of freedom, have been combined in order 
to construct antisymmetric basis functions. In fact, the HH functions, as normally defined, 
do not have well defined properties under particle permutation; this results from the selection 
of a particular ordering of the particles in the definition of the Jacobi coordinates and, as 
a consequence, of the hyperangular coordinates used to define the HH functions. Changing 
the ordering of the particles, it is possible to define a new set of Jacobi coordinates and, 
accordingly, HH functions depending on the hyperangular variables obtained from this new 
set. To be noticed that the HH functions defined using a particular choice of the Jacobi 
coordinates form a complete basis. 

The HH functions defined in one set of Jacobi coordinates can be transformed to HH 
functions defined in another set. In this transformation (permutation) the grand angular 
quantum number K, which identifies a subset of HH functions, is conserved. For finite 
values of K, the dimension Nk of this subset is finite, and therefore a finite number of 
HH functions, having all the same value of K, are necessary to describe a HH function 
having the same value of K but defined in a different Jacobi set. The coefficients of the 
transformation can be collected in a matrix having the dimension Nx for each number of 
particles. For A = 3 these matrix elements are the Raynal-Revai coefficients [3]. For A > 3 
the coefficients cannot be given in a close form, and a few methods have been devised for 

the, c alcula t lons aaaj The MAge * ^ for the ^ 

of basis elements with well defined permutational symmetry. In fact, each subset defined 
by K is invariant under particle permutation, as a consequence, the constructions of basis 
elements with that property is performed as linear combinations of HH functions having the 
same value of K. Different schemes to construct hyperspherical functions with an arbitrary 
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permutational symmetry are given in Refs. |6|, |8|, |St lldl . 111] . Recently, a procedure for 



constructing HH functions in terms of a single particle basis has been proposed in Ref. [12]. 

In problems in which the A-body system is composed by identical particles the wave func- 
tion of the system has to be completely symmetric or antisymmetric in the case of bosons 
or fermions, respectively. Considering a Hilbert space extended to spin and isospin degrees 
of freedom, the construction of HH functions having well defined permutational properties 
allows for a reduction of the large degeneracy of the basis. In general the completely sym- 
metric or antisymmetric basis functions are a small part of the total Hilbert space. However, 
the difficulties of constructing HH functions with well defined permutational symmetries in- 
crease with A and K; therefore, the preliminary step of constructing basis functions with 
well defined permutational symmetry for A particles could be sometimes very difficult to 
carry out. 

In the present paper we investigate a different strategy. We intend to perform the de- 
scription of a A-body system using the HH basis defined on one set of Jacobi coordinates, 
the reference set, and not having a well defined behaviour under particle permutation. We 
will loose the advantage of using a reduced part of the total Hilbert space; however, we will 
gain in simplicity in the calculation of the matrix elements. By including all HH basis ele- 
ments up to a certain grand angular momentum K, the diagonalization of the Hamiltonian 
matrix will produce eigenvectors reflecting its symmetries. If the Hamiltonian commutes 
with the group of permutations of A objects, Sa, in the case of non degenerated eigenvalues, 
the eigenvectors will have a well defined permutation symmetry, and can be organized in 
accordance with the irreducible representations of Sa- Therefore, identifying those eigen- 
vectors with the desired symmetry, the corresponding energies can be considered variational 
estimates. In particular, it will be possible to identify a subset of eigenvectors and eigen- 
values corresponding exactly to those that would be obtained performing the preliminary 
symmetrization of the states. The disadvantage of this method results in the large dimension 
of the matrices to be diagonalized. However, at present, different techniques are available 
to treat (at least partially) this problem. 

For a system interacting through a two-body potential V(i,j), the potential energy op- 
erator results in a sum over pairs. Its matrix elements can be reduced to one term, let us 
say V(l, 2), times the number of pairs when symmetric or antisymmetric state functions are 
considered. When HH functions without well defined permutation behaviour are used, the 
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calculation of the potential energy operator cannot be reduced to the computation of one 
term. So we have to face the problem of computing the matrix elements of a general term 
V(i,j) between HH functions defined in the reference set of Jacobi coordinates in which the 
distance between particles has not a simple form. 

The calculation of V(i,j) in the reference set of Jacobi coordinates is performed in two 
steps: (i) first, we use a property of the HH basis which allows to expand a general function 
of the coordinates in terms of a subset of the basis called the potential basis (PB) 
(ii) then, as for the case of a generic HH basis element, each PB element is transformed 
to the HH basis defined in the reference set of Jacobi coordinates. In the case of the PB, 
the transformation coefficients are known analytically for each value of K and for a general 
number of particles A. In this way, each term V(i,j) of the potential energy can be first 
expanded in the PB, and then transformed to HH functions defined in the reference set. So, 
after this procedure, the potential energy will be expressed in terms of HH functions. As we 
will see, the computation of the matrix elements of the potential energy is now very simple 
since it results in a combination of integrals of three HH functions. A compact form suitable 
for a numeric treatment of the problem will be given. 

The derivation and implementation of the final expression for the potential energy in the 
calculation of bound states is the main subject of the present paper. As a simple application, 
a system of three and four nucleons interacting through a central potential will be analyzed. 
Different symmetries will appear considering or not the Coulomb interaction between two 
protons. To be noticed that when antisymmetrized basis functions are used to describe 
three or four nucleons, the presence of the Coulomb interaction implies that states with 
total isospin T = 1/2,3/2 (for A=3) and T = 0,1,2 (for A = 4) have to be considered. 
The extension of the Hilbert space to include these terms increases the dimension of the 
problem resulting comparable to that one in which the antisymmetrization of the basis is 
not performed. Finally we would like to stress that the present paper is the first step in a 
program devoted to applications of the HH basis to systems with A > 4 interacting through 
realistic potentials. 

The paper is organized as follows, section [Til is devoted to a brief description of the HH 
basis. In sections lIHI the expression for the potential energy in terms of HH states are given. 
In section IIVI the results for the examples proposed are shown. Section [V] includes a brief 
discussion of the results and the perspectives of the present work. 
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II. THE HARMONIC HYPERSPHERICAL BASIS 



In this section we present a brief overview of the properties of the HH basis following 
Ref. We start with the following definition of the Jacobi coordinates for an A body 
system: 



where m is a reference mass, N = A — 1 and we have defined 

j ^ j 

Mj = J2 m i' X i = Jf.J2 miVi - ( 2 ) 

i=l 3 i=l 

Let us note that if all the masses are equal, rrii = m , Eq. ([!]) simplifies to 



x w- 3 +i = \ I -^p-y fo+i - Xj) , j = 1, . . . , N . (3) 
For a given set of Jacobi coordinates x 1; . . . , x^r, we can introduce the hyperradius p 



o=(£*0 = 2 £<"■•- x > 2 = iEfe-^.) 2 • («) 

and the hyperangular coordinates Qjv 

OjV = (#1, • • -,XN,<j>2, ■ ■ ■ >0JV) , (5) 

with the hyperangles <pi defined via 

cos d>i — — % , i — 2,...,N. (6) 

^xj + --- + xf K) 

The radial components of the Jacobi coordinates can be expressed in terms of the hyper- 
spherical coordinates 

Xn = P COS (j) AT 

xn-i = psin^Ar cos0at_i 

(7) 

Xi = p sin • • • sin (f) i+ i cos 
X\ = p sin (f)pf • • • sin 3 sin </> 2 . 
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Using the above hyperspherical angles f^, the surface element becomes 



N 



dfl^ = sin 61 dOi dipi sin 9j d9j dipj (cos 4>j) 2 (sin (f)j) 3 -' 4 d(f>j , 

J'=2 



(8) 



and the Laplacian operator 



2 , a 2 3jv-i a A^(!i K ; 



i=l 



<9p 2 p <9p p 2 



(9) 



where the A 2 v (fij V ) is the generalization of the angular momentum and is called grand angular 
operator. 

The HH functions y[K](&N) are the eigenvectors of the grand angular momentum operator 

A 2 N (n N ) + k{k + 3N - 2)) y [K] (n N ) = o . (io) 

They can be expressed in terms of the usual harmonic functions Yi m {x) and of the Jacobi 
polynomials P®' b (z). In fact, the explicit expression for the HH functions is 



N 



' N 



N 



J=2 



(11) 



where [X] stands for the set of quantum numbers l±, . . . , In, mi, . . . , m^, n-i, . . . , tin, the 
hyperspherical polynomial is 



(12) 



where the K~ quantum numbers are defined as 



K j = Y^{l i + 2n i ), ni = 0, K = K N , 



(13) 



i=i 



and the normalization factor 



2i/ J -r(i/ j -nj)nj! 



r(^. _ „ . _ i . _ 1/2) r( nj + ij + 3/2) ' 



(14) 



with v,j = Kj+3j /2 — 1. The quantum number K is also known as grand angular momentum. 
The HH functions are normalized 



(15) 
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moreover, the HH basis is complete 



[K] V ' 



5 



3JV-l/o' 



(n 



A' 



AT, 



(16) 



With the above definitions, the HH functions do not have well defined total orbital angular 
momentum L and z projection M. It is possible to construct HH functions having well 



defined values of LM by coupling the functions Y\ 



ij-rrij y x jj 



This can be achieved using 



different coupling schemes. Accordingly we can define the following HH function 



[K] 



A J 



J LM 



A 



J'=2 



(17) 



having well defined values of LM, although the particular coupling scheme is not indicated. 
The set of quantum numbers [K] now includes the values of LM plus N — 2 intermediate 
/-values instead of the iV magnetic numbers rrij. When necessary, the explicit coupling 
scheme of the above HH function will be given. 



A. Potential Basis 



If we have a function which depends only on the difference of two particle positions, 
/(ri — r 2 ), we can use a subset of the HH's to expand that function, called the potential 
basis (PB) [l|. Let's introduce the Jacobi coordinates such that = ri — then the PB 
subset is defined by the following condition (f2i 2 = (xm,<Pn) and Qn = (Qn-i,xn,4>n)) 

A5r_ 1 (n^-i)7&(n 13 ) = J (18) 

where (n,l,m) = (njv, In, itin), an d by 

^U^V^u) = —K(K + 3N- 2)V^ +l (n 12 ) , (19) 

with K = I + 2n. Thus, the PB is a subset of the HH's which depends only on (xjv,0jv) 
variables, and which is specified by only three quantum numbers n,l,m, instead of the 
3N — 1. The PB basis element has well defined angular momentum I and projection m. The 
expression of the PB elements is: 

V l ^ +l (n 12 ) = Y m (x^)(cos^)' J R5 Ar - 1 )/ 2 - 1 ' z+1 / 2 (cos20 Ar )r [o] ( J D - 3) , (20) 
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where (D = 3JV) 



Y W (D-3) 



r[(D-3)/2] 

2 7r (D-3)/2 



1/2 



is the normalization verifying 



[y [0] (£>-3)] 2 dn 



AT-l 



1. 



The surface element is conveniently written as 



(21) 



(22) 



(23) 



We can extend the definition of the PB elements to depend on the coordinates of a 



general pair as ^2n+i(^y)- The coordinates Qij = (xn,(Pn) are now defined by a 

different ordering of the particles entering in the Jacobi coordinates such that xjv = r« — rj. 
One important property of the PB elements is the following. When a PB element is defined 
in the space spanned by the coordinates Q^, its expression in terms of HH functions defined 
in the reference set corresponding to the ordering of the particles 1,2, ....,N is known 
and results 



[K'=2n+l] 

where the coefficients ^C^^{<f^) are given by the following relation 

The angles (p^ = {p 1 ^, ■ ■ ■ , p>i} defined from the following kinematic rotation vector 

z((p^) = xat cos </5^+xjv-i sin (p^ cos <p%_x + ■ ■ • + xi sin tp 1 ^ sin f^^-i • • ■ sin ^ (26) 

are chosen to verify z(ip 1 ^) = rj — rj. The hyperangles Q l J are defined as Q l J = 
{x, ... x, ipx, . . . , p> % i}, with x repeated iV-times. The particular form of the HH function 
y^nipfi) produces the coefficients of Eq. fT25l) to be independent of m. In Eq. fT24|) the sum 
over all quantum numbers [K'\ is limited by the condition 2n + 1 = K', showing that a PB 
basis element depending on fl^ can be given as a linear combination of HH basis elements 
having the same value of grand angular quantum number but depending on Q^. A generic 
function /(r^ — r^) can be expanded in terms of the PB as 

f{Ti - Tj) = E fnlmiP^n+l^ij) , (27) 
nlm 
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with 



fnlm(p) = J dVtijf^i-Tj) J dQ N -l[V^ +l {%)Y 



Y (D 



J dSlqlfa (cos^) / P n 3 ^- 1 )/ 2 - 1 '' +1 / 2 (cos2^)/(r 4 - r,) . 



(28) 



The functions f n im{p) are the hyperradial multipoles. Introducing the transformation of 
Eq. (l24|) in the above expressions it is possible to write a general function /(r, — r^) in terms 
of HH functions given in the reference set. We will use this property for the potential energy 
of an A-body system. 

III. THE POTENTIAL ENERGY IN TERMS OF HH FUNCTIONS 

A local two-body interaction can be put in the form 



V,(r. 



where we use the compact notation 



A h {i,j)®Y h {h 



LM 



^2 (lim 1 l 2 m 2 \LM)A hmi (i,j)Yi 2m2 (f ij ) . 



(29) 



(30) 



m\m,2 



Aim(i,j) is an operator independent of the coordinates r^, and the coupling with the spher- 
ical harmonics to zero in Eq. (l29|) shows that the potential is a scalar in total space. We can 
use the PB elements to expand each /-term of the expansion 



V(i,j) = £ [Mi J) ® 'PLhW] VI{r) , 



(31) 



In 



where the functions V^(p) are obtained from the following integral in the hyperangular space 

•)l,m 

(32) 



V&>) = J d^V^YUhi) I dQ^V^Q^r 
1 



Y (D - 3) 
The complete potential energy is 



^(cos^) 2+ '(sin^) 3JV - 4 P n 3 ^- 1 )/ 2 - 1 '' +1 / 2 (cos2^)^(r lJ ) 



J2V(i,j) = EE iMhj) ®K + i(^)] v n(p) ■ 

i<j i<j In 



(33) 



It would be convenient to have the potential energy expressed in the coordinates defined 
by Q (in the following we drop the suffix N for the reference set). To this end we transform 
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the PB elements obtaining 



i<j In [K'=2n+l] i<j 

= E F »^), 

In 



where we have defined 



= E (E Wc &(^)) [^i)®^](o)] . (35) 

The final form of Eq. (|3^|) gives a general expression for the potential energy in terms of 
the HH basis elements. In the case of central potentials I = and A\ m = 1, and the above 
expressions reduce to (omitting the indices Z = 0,m = 0) 

E y (^') = EE( E {N) c^ ij )yiK'm)v n ( P ) 

i<j i<j n [K'=2n] ' 

= Ew E (E (iY) ^'](^'))^](fi) (36) 

n [K'=2n] ^ i<j ' 



with 



n 

= E (E%(^ 

[Jf'=2n] ^ i<j ' 
= E^n(%) • 



(n) 



(37) 



The matrix elements of the potential energy between two different HH basis elements 
result 

i<j nl 

The above expression represents an integral in the hyperangular space, and shows the tensor 
product form between the hyperradius and the hyperangular coordinates which is typical 
using the HH basis. The matrix elements of the operators G l n {Q) are independent of the 
potential 



(y^mmy$T)«= E E w ^> iJ ) 



[K'=2n+l] %<j 



X 



E 7m=f Ai - m ^ 3) J dn y l m m yffi m ■ 
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Each Q n is a combination of HH functions with grand orbital momentum K' = 2n + I, 
therefore its matrix elements follow a triangular relation. In fact, given K\ and K 2 , the 
values of n, I to be considered in the sum of Eq. (138|) are limited by the relation \K\ — K 2 \ < 
2n + I < Ki + K 2 . A triangular relation is also verified by the orbital angular momenta: 
\Li — L 2 \ < I < L\ + L 2 . Furthermore, the matrix elements of Q n includes the computation 
of integrals of three HH basis elements: 



These integrals factorize in products of one-dimensional integrals consisting of either three 
hyperspherical polynomials or three spherical harmonics that can be obtained analytically 
or very efficiently using quadratures. 

As shown in Eq.( l35|) . each function G n (Q) is symmetric in the particle indices, therefore its 
corresponding eigenvectors will have well defined symmetry under particle permutations. For 
example, when A lm (i,j) = 1, K x = K 2 = K, I = 0, implying (L u M 2 ) = (L 2 , M 2 ) = (L, M), 
and 2n = 2K, the following elements 



form a matrix by varying all the quantum numbers in [K] with fixed values of K and L, M. 
The dimension of the matrix is given by all HH functions with grand angular quantum 
number K coupled to L, M. Its eigenvectors, which are combinations of this family of HH 
functions, will have well defined permutational symmetry. This reflects the fact that each K- 
subset is invariant under particle permutations. Therefore, the diagonalization of the above 
matrix is a way to construct basis elements with well defined permutational symmetry using 
HH functions with fix values of K and L. 

IV. APPLICATION TO SYSTEMS WITH A = 3, 4 

In the description of bound states in an A-body system it is common to use basis elements 
having the required symmetry, symmetric states for bosons or antisymmetric states for 
fermions. In the present section we will analyze the use of the HH basis without the initial 
symmetrization or antisymmetrization of the basis. Although the basis elements have not the 
required symmetry, the eigenvectors of the Hamiltonian will have a well defined symmetry 




(40) 




(41) 
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reflecting the symmetries appearing in the Hamiltonian. Therefore, among all eigenvectors 
and eigenvalues, the physical states have to be identified. 

By taking opportune linear combinations of the HH basis elements, specific symmetries 
under particle permutation can be constructed for fixed values of K. Therefore two calcu- 
lations, one in which all HH states up to a maximum value of K are considered and the 
other in which states with a particular symmetry up to the same value of K are considered, 
produce the same eigenvectors and eigenvalues. Of course, in the first calculation will appear 
eigenvectors and eigenvalues belonging to other symmetries not present in the second cal- 
culation. The simplification of avoiding the initial basis symmetrization is counterbalanced 
by the larger dimension of the Hamiltonian matrix. 

Limiting the discussion to central potentials, Eq. (1381) is well suited for a direct application 
of the HH basis. Let us introduced the following orthonormal basis element 



( P Q\m[K})= (/3 (a+1)/2 



ml 



(a + m) 



(42) 



where Lm\/3p) is a Laguerre polynomial with a = 3N — 1 and (3 a variational non-linear 
parameter. We will discuss the case L = for A = 3,4. The HH basis elements are 



Yi{xi) <g>Yi(x 2 ) 



(43) 



for A = 3, and 



\Yi 1 (x 1 )®Y la (x2)] h ®Y h (x 3 ) 



(44) 



J o 



for A = 4. The corresponding matrix elements of the Hamiltonian are 



(to' [K']\H\m[K]) = (m'[K']\T + V\m[K]) 



(45) 



The matrix elements of the potential energy corresponding to the fi-space have been dis- 
cussed in the previous section. Integrating also on p-space they result 

(m'[K']\V\m[K}) = ^(m'|K(p)| m) p { [K']\g n (n)\ [K\) a = £TOJ(0fa [it] ) ■ (46) 
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The matrix elements of the kinetic energy are the following 



h 2 N 



T K ' m ';Km = (rri [K'\ I -_VV^|m[K]) 



m 
i=i 



— {m'[K'}\ — + — + -JDU.\ m [K\) 

m dp 2 p dp p 2 (47) 

— °[K\\K'\J- m >, m 



— [T« m - K{K + 3N - 2)T® 



with 



- L rn'rn , iJ m.m / I 



1 „ / m'! / ml 



i ' V (a + m')! V (ot + m)\ 



00 



'm,0; 



(48) 



and 



r - ^ \IW^.]]T^. [ ^e-^L £ > W (l)L£» W . (49) 

Using the properties of the Laguerre polynomials, these integrals can be calculated analyti- 
cally. 

Therefore, the matrix elements of the Hamiltonian are sums of tensor products of two 
matrices, one calculated on p-space, depending on indices m, m', and one calculated on 
fi-space, depending on the indices [K], [K'\ 

(m'[K'\\H\m[K\) = —^—^(T^} m —K . (50) 

n 

If we introduce the diagonal matrix D such that ([K 1 ] \ D \ [K]) = 5[kuk'\K{K + 3iV — 2), 
and the identity matrix / in .fT-space, we can rewrite the Hamiltonian schematically as 

H = -^-(1 ® WT + D ® ^T) +YQ n ®V n , (51) 

in which the tensor product character of the expression is explicitly given. A scheme to 
diagonalize such a matrix is given in the Appendix [A] 

In the following we give results for nucleon systems with A = 3, 4 using the Volkov 
potential 

V(r) = V R e~ r2/R i + V A e~ r2/R2 i (52) 
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with V R = 144.86 MeV, R x = 0.82 fm, V A = -83.34 MeV, and R 2 = 1.6 fin. The nucleons 
are considered to have the same mass chosen to be equal to the reference mass m and 
corresponding to h 2 /m = 41.47 MeVfm -2 . With this parametrization of the potential, the 
two-nucleon system has a binding energy of 0.54592 MeV. 

This potential has been used several times in the literature making its use very useful 
to compare different methods [ll, Q, Q]- The results will be obtained after a direct 



diagonalization of the Hamiltonian matrix of Eq. fl50|) including m max + 1 Laguerre poly- 
nomials with a fix value of /3, and all HH states corresponding to maximum value of the 
grand angular momentum K max . The scale parameter (3 can be used as a non-linear param- 
eter to study the convergence in the index m = 0, 1, . . . ,m max , with m max the maximum 
value considered. In the present analysis the convergence will be studied with respect to 
the index K max , therefore, the number of Laguerre polynomials at each step, m max + 1, 
will be sufficiently large to guarantee independence from (3 of the physical eigenvalues and 
eigenvectors. 

In Table U we show the different symmetries of the eigenvectors and the corresponding 
eigenvalues, for A = 4, in the particular case in which the Hamiltonian matrix has been 
diagonalized for m max — 0, f3 — 2 fm -1 , and K max = 6. In this case the total dimension of 
the matrix is 56 with 32 "even" elements, corresponding to even values of 1%, and 24 "odd" 
elements corresponding to odd values of Z3. In particular there are 6 totally symmetric states, 
irreducible representation [4] using the Yamagouchi symbol, 2 totally antisymmetric states, 
[l 4 ], 8 states belonging to the three-dimensional irreducible representation [31], 6 states be- 
longing to the two-dimensional irreducible representation [2 2 ] , and 4 states belonging to the 
three-dimensional irreducible representation [2 l 2 ] . The lowest eigenvalue of each irreducible 
representation is given in the table. In the last two columns of the table, the eigenvalues are 
reported considering separately even and odd basis elements. Symmetric states are formed 
exclusively by even-basis element whereas antisymmetric states are formed exclusively by 
odd-basis elements. The three mixed symmetries, one two-dimensional and the other two 
three-dimensional, show degenerate eigenvalues. In order to distinguish between the two 
three-dimensional mixed symmetries, we observe that the three degenerate eigenvalues di- 
vide differently in even and odd elements. The mixed symmetry [3 1] is twice degenerate 
when the expansion basis is restricted to even states whereas the mixed symmetry [2 l 2 ] is 
not. Therefore performing two different diagonalizations, one using a restricted basis con- 
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sidering only even states and one considering only odd states, all the symmetries can be 
identified. 

Furthermore, we can see from Table [I] that a bound state appears in correspondence to a 
symmetric state. The fact that only one spatial symmetry is present in the bound state is a 
direct consequence of using a central potential. The final antisymmetrization of the state, as 
required in the case of four nucleons, is performed by multiplying the spatial symmetric wave 
function by the corresponding spin functions, singlet spin states S12 = for the two protons 
labelled (1,2) and S34 = for the two neutrons labelled (3,4). In the case of using the 
isospin formalism, the spatial symmetric state is multiplied by a four nucleon antisymmetric 
spin-isospin function having total spin S = and total isospin T = 0. 

In Tables ITU and [TTT1 the convergence of the ground state binding energies for A = 3, 4 are 
given as a function of K max , respectively. In the last column the point Coulomb interaction 
between the two protons, labelled as particles (1, 2), has been considered. In the case without 
the Coulomb potential, spatial component of the ground state is completely symmetric. 
When the Coulomb potential is taken into account this component is symmetric with respect 
to particles (1,2). For A = 4, it is also symmetric with respect to the particles (3,4), the 
two neutrons. In this case it is convenient to introduce the if-type Jacobi coordinates (for 
a recent application see Ref. Q): 

x 3 = r 2 - n 

_ r 4 + r 3 r 2 + n 

X2 " ~7I 7T (53) 

xi = r 4 - r 3 , 

and construct HH basis elements based on this type of coordinates. These HH functions 
are linear combinations of the HH function based on the K-type coordinates, introduced 
in Eq.Q and used in the previous sections, at fixed values of the grand angular quantum 
number K. As example, in Table HTT1 the two different types of Jacobi coordinates have been 
considered. The dimension of the bases indicated corresponds to taking into account even 
basis elements which are the only ones entering in the construction of the bound states. As 
stated before, for the if-type Jacobi coordinates this means to take even values of I3. For 
the if-type, both l\ and I3 are taken even. The dimension of the problem for obtaining the 
eigenvalue at K max = 30 results to be 72 for A = 3, and 7872 (4056) for A = 4, using the 
fT-type (if -type). The use of the if -type Jacobi coordinates reduces the dimension of the 
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problem by nearly a factor of two. 

The calculations corresponding to the two different types of coordinates differ in the set of 
angles ip^ defined in Eq. fl26j) reflecting the different way of defining the interparticle distances 
in both cases. In the case in which the symmetric states are identified and constructed 
before diagonalization, the dimension is reduced to 27 for A = 3 and around 600 for A = 4. 
We observe a considerable reduction in the dimension of the eigenvalue problem for the 
symmetrized basis. However the computational cost of constructing HH states with specific 
permutational symmetry has to be compared to the simplicity in constructing the matrix 



elements of Q n and in solving t 
symmetrized HH basis of Ref 



re system of Eq.(l5"Tl). To be noticed that the results using the 



coincide with the results presented here for each value 



of K mnr . For the sake of comparison, the results using the stochastic variational method 



(SVM) 13J are shown in the table. 

When the Coulomb potential between protons is included the system can be treated as 
composed by two different species, the protons and the neutrons, having different interactions 
and slightly different masses. Using the complete HH basis this cause no extra difficulties 
since the following term can be added to the Hamiltonian 

(54) 

n 

where V^ m are the hyperradial matrix elements of the Coulomb potential multipoles and 
3~Jk']\K] * s a ma tri x equivalent to ^[^m^i with the only difference that the sum over (i, j) is 
limited to protons. This term has the tensor product form and therefore the Hamiltonian 
reads: 

H = I® (1) T + D ® (2) T + <3n ® V n + ® V° . (55) 

n n 

In the above equation protons and neutrons are assumed to interact with the same short 
range potential. For realistic potentials this is not the case and this assumption can be re- 
laxed dividing the potential energy in three parts, one for the interaction between protons, 
one for the interaction between neutrons and one for the interaction between protons and 
neutrons. To be noticed that using the complete HH basis the dimension of the problem does 
not change by distinguishing protons and neutrons or not. However, as is a common proce- 
dure, it is possible to treat the system as composed by identical particles using the isospin 
formalism. The Coulomb potential breaks the isospin symmetry and the use of antisymmet- 
ric states requires the inclusion of different isospin components in the wave function. The 
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A = 3 bound state will have isospin T = 1/2, 3/2 components whereas the A = 4 bound state 
will have T = 0, 1, 2 components. After including all these components the two procedures, 
one using the complete HH basis and the other using antisymmetrized states, will produce 
the same eigenvalues. An example for this case is given in Table IV in which the results for 
A = 3,4 using antisymmetric basis states, including the different isospin components, are 
shown. For A = 3, the T = 1/2 component is by far the most important one; however, the 
exact result is obtained after including both components, T = 1/2 and 3/2. To be noticed 
that the dimension of the basis using antisymmetrized HH states with isospin components 
T — 1/2, 3/2 is the same of the complete HH states using even basis elements. Therefore in 
this case the preliminary antisymmetrization of the basis is not convenient. For A = 4 the 
T = component is by far the most important; however, the exact result is obtained after 
including the three isospin components T = 0, 1,2. In this case the dimension of the basis 
using even HH states up to K max is greater than that using antisymmetrized basis elements, 
since in the symmetrization with respect the two neutrons is not included automatically in 
the even HH states and has to be constructed by the diagonalization procedure. However 
the difference in the dimension of the two cases is considerably reduced with respect to the 
case in which the Coulomb potential was not included. 

The equivalence of the last columns of Tables II and III with columns third and sixth of 
Table IV illustrates the simplicity of treating symmetry breaking terms using the HH basis 
without permutational symmetry. 

V. CONCLUSIONS 

In this work we have presented a direct use of the HH basis in the description of a A- 
body system. The basis has neither been symmetrized nor antisymmetrized as required by a 
system of identical particles. However, the eigenvectors of the Hamiltonian have well defined 
permutation symmetry. Among all the eigenvectors, the physical ones can be identified. The 
benefit of the direct use of the HH basis is based on a particular simple form used to represent 
the potential energy. Each term of the two-body potential V(i,j) has been expanded in the 
potential basis and then expressed in terms of the HH basis, defined in the reference set, by 
using the corresponding transformation coefficients. These coefficients are known for each 
value of K and for a general number of particles A. Once the potential has been expressed 
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in terms of the HH basis, it results in a sum of tensor product terms originated from the 
separation of the hyperradial and the hyperangular coordinates inherent to the method. 
Moreover, the kinetic energy can be put in a tensor product form too. Therefore, the matrix 
representation of the Hamiltonian is expressed as a sum of tensor product matrices, and 
this particular form can be diagonalized very efficiently using the technique given in the 
Appendix. As a test case we have studied three and four nucleons interacting through a 
central potential, the Volkov potential, used many times in the literature. We have shown 
how the symmetries are present in the spectrum and can be identified. The symmetric and 
antisymmetric states appear as singlets, whereas the mixed symmetries appear as multiplets. 
We have identified all symmetries by dividing the spectrum in even and odd components. 
In the studied cases, only one bound state appears for A = 3 and 4 corresponding to a 
symmetric state. To be noticed that if the potential depends on the spin- isospin degrees 
of freedom, the Hamiltonian will still present the tensor product form in the hyperradial, 
hyperangular, spin and isospin spaces. 

For A = 4 we have solved the problem using two different types of Jacobi coordinates, 
namely the K-type corresponding to a 3+1 configuration, and the if -type corresponding 
to a 2+2 configuration. The calculations using one or the other set differ in the values of 
the angles (p 1 *, which can be considered as input parameters. Therefore the method gives a 
systematic way of introducing the different types of Jacobi coordinates. The convenience of 
selecting one specific type is related to its capability to produce basis states having partially 
the required symmetry with a reduction of the total dimension of the problem. In the cases 
presented here, the A = 4 bound state is constructed using basis elements based on the 
fT-type Jacobi coordinates with even values of 1% or based on the if -type with both li and 
1% restricted to even values. The latter resulting in a basis with a dimension smaller by a 
factor of two. 

A further benefit of using the complete HH basis is obtained when symmetry breaking 
terms are included in the Hamiltonian. The complete basis will generate eigenvectors having 
specific permutation symmetries reflecting the symmetries present in the Hamiltonian. The 
complexity of the numerical problem does not increase when these terms are present. This is 
not the case when symmetrized or antisymmetrized basis are used. For example, in the case 
of a nuclear system, the presence of charge symmetry breaking terms requires the extension 
of the basis to include all the isospin components. As a specific example, here we have 
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analyzed the case of the Coulomb interaction between protons. The results using the HH 
basis without well defined permutation symmetry have been compared to the case in which 
antisymmetrized HH basis have been used. In the latter case the different isospin components 
entering in the wave function have to be included, resulting in spatial components having 
more than one symmetry. Accordingly the dimension of the basis increases. To this respect, 
the numerical effort to reduce the Hilbert space to subspaces with specific permutation 
symmetry is discussed in Ref.js]. When one spatial symmetry is required, as for example a 
completely symmetric spatial state, the convenience of constructing symmetric HH state is 
obvious. When several spatial symmetries are present in the wave function, as in the case of 
an A-nucleus wave function, the convenience of constructing HH states with different spatial 
symmetries has to be compared to the capability of solving a large eigenvalue problem, for 
example that one given in Eq. (l51D . 

nnri 

In Refs. [14|, [17|, [18| the HH basis, used to describe three- and four-nucleon bound states, 



is antisymmetrized in the following way. The total wave function is expanded in angular- 
spin-isospin channels and, for each channel, it is written as a sum of Faddeev-like amplitudes, 
each of them antisymmetric in the pair (i, j). In this way the total wave function results 
antisymmetric. Then, each -amplitudes is expanded in the HH basis defined from 
Jacobi vectors corresponding to the different ordering of the particles. As a consequence, the 
amplitudes for the different channels are not orthogonal, resulting in a non-orthogonal basis. 
For large values of K the non orthogonality of the basis could causes numerical instabilities. 
In particular, for A = 4, this problem is overcome performing an orthonormalization of 
the basis using the Gram-Schmidt technique with quadruple precision in the numerical 
treatment of the process. Therefore, the extension to A > 4 systems appears to be difficult. 
On the other hand, the direct use of the HH basis without antisymmetrization circumvents 
this problem. Therefore, the method presented here has to be considered a first step in 
a program devoted to the application of the HH basis to systems with A > 4. Further 
works along this line are the extension of the method to treat realistic interactions and the 
numerical implementation of the Hamiltonian of Eq.( l5"Tj) to systems with A = 5,6. The 
extension of the method to treat three- nucleon interaction terms is also possible. In fact, 
the transformation of the spatial part of a three-nucleon interaction W(i, j, k) in terms of HH 
functions constructed in the reference set can be performed using the algorithm developed 
for the "triplet basis" in Ref. j^. 
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APPENDIX A: EFFICIENT MATRIX- VECTOR PRODUCT FOR TENSOR- 
PRODUCT MATRICES 



The algorithm we used to diagonalize the Hamiltonian is an iterative one, namely Lane- 

n 

zos's [19| . These kind of algorithms are useful whenever an efficient matrix-vector product 
can be used, as in the case of sparse matrices; in the specific calculation, we have the product 
between a tensor-product matrix M = A\ ® A 2 , and a vector v 

w = M ■ v = (A 1 <g> A 2 ) ■ v, (Al) 

with Ai a n x n matrix, A2 a m x m matrix, and v a (n ■ m) -dimensional vector. 

The product is done in three steps: (i) first, the vector v is reshaped inamxn matrix 
V; (ii) then, the following matrix products are performed 

W = {Ai- (A 2 ■ V) T ) T ■ (A2) 

(iii) finally, the matrix W is reshaped into the (n ■ m)-dimensional vector w, which is the 
result of the multiplication. The above algorithm is easily generalized to tensor-products of 
A;-matrices 20]. 
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TABLE I: Lowest Volkov-energy eigenvalues of each irreducible representations of S4 for the N = 4 
case, with m max = 0, K max = 6, and [3 = 2fm _1 . The multiplets are further identified as being 
symmetric or anti-symmetric under permutation of particles 1-2. 



Irreps 


Eigen's (MeV) 


Sym(l-2) 


AntiSym(l-2) 


[4] 


1 1 1 2 1 3 1 4 




-25.794 


-25.794 




[2 2 1 


1 

3 


2 
4 


27.680 
27.680 


27.680 


27.680 


[3 1] 


1 

4 


2 1 3 




28.430 
28.430 
28.430 


28.430 
28.430 


28.430 


[21 2 ] 


a 
4 




102.85 
102.85 
102.85 


102.85 


102.85 
102.85 












[I 4 ] 




L 

> 
i 
1 


199.56 




199.56 
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TABLE II: Results 


for the Volkov's potential, 


as a function of K max using 


30 Laguerre's polynomi- 


als, and /3 = 3 fm 1 for the three-body case. 


In the last column the result 


s including 


the Coulomb 


potential are given. 










K 


/Vtttt 

iv HH 




R CMpVI 

J-J I 1VJ.C V J 







1 


7.7075 




fi QQ9fi 


10 


12 


8.4157 




7.7083 


20 


36 


8.4623 




7.7566 


30 


72 


8.4647 




7.7693 


40 


121 


8.4649 




7.7694 



SVM[13] 30 8.46 



TABLE III: A = 4 results for the Volkov's potential, using 25 Laguerre's polynomials, and j3 = 
2fm _1 . Two different types of Jacobi coordinates have been used. In the last two columns the 
results without and with Coulomb potential, using independently K- or //-type Jacobi coordinates, 
are given, respectively. At fixed value of K max the results, using either one or the other type of 
coordinates, coincide. 

ifmax iVHH(if-type) N UH {H-type) E (MeV) 

1 1 28.580 27.748 

10 136 78 30.278 29.456 

20 1547 819 30.416 29.596 

30 7872 4056 30.420 29.599 

SVM[13] 50 30.42 
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TABLE IV: Contributions to the bound state energies, for A = 3,4, of the different isospin com- 
ponents using antisymmetrized HH functions 



A = 3 A = 4 



K 


T = 1/2 


T = 1/2,3/2 


K 


T = 


T = 0,1,2 





6.9926 


6.9926 





27.748 


27.748 


10 


7.7072 


7.7083 


10 


29.453 


29.456 


20 


7.7555 


7.7566 


20 


29.594 


29.596 


30 


7.7582 


7.7593 


30 


29.596 


29.599 


40 


7.7583 


7.7594 
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